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Comment  on  “A  Coupled  Dynamic-Thermodynamic  Model  of  an  Ice-Ocean 
System  in  the  Marginal  Ice  Zone”  by  Sirpa  Hakkinen 

P.  C.  Chu  and  Roland  W.  Garwood,  Jr. 

Department  of  Oceanography,  Naval  Postgraduate  School,  Monterey,  California 


The  recent  article  by  Hakkinen  [1987]  creatively  addresses  a 
number  of  important  dynamic  and  thermodynamic  processes 
in  the  marginal  ice  zone  (MIZ).  We  believe  that  her  study  will 
stimulate  further  research  on  the  coupled  ice-ocean  system. 
For  now,  we  would  like  to  offer  the  following  discussion  of  the 
need  to  augment  the  entrainment  equation  to  include  surface 
buoyancy  flux  for  general  application  to  the  MIZ. 

Inclusion  of  Buoyancy  Flux  in  Entrainment 
Formulation 

In  Hakkinen’s  study  of  the  coupled  ice-ocean  system  the 
ocean  is  considered  to  have  two  layers.  An  essential  part  of  the 
thermodynamics  is  the  mixed  layer  entrainment  of  the  lower 
layer  by  the  surface  layer.  However,  entrainment  was  parame¬ 
terized  using  only  the  friction  velocity  at  the  ocean  surface, 
including  an  ad  hoc  e-folding  depth  dependence  from  Pollard 
et  at.  [1983]: 

=  6 u„3  exp  ( — Ai,//j0)/(c„2  +  g*hj  (1) 

where  fi,  is  the  upper  layer  depth,  is  the  water  surface 
friction  velocity,  g*  is  the  reduced  gravity,  h0  is  the  length 
scale  for  e-folding  damping,  and  Cm2  is  a  parameterization  for 
unsteadiness  in  the  turbulent  kinetic  energy  (TKE)  budget. 
The  values  of  h0  and  Cm  are  specified  by  Hakkinen  to  be  20  m 
and  0.03  m/s,  respectively.  Equation  (1)  is  insufficient  for  gen¬ 
eral  application  to  the  MIZ  because  it  does  not  include  buoy¬ 
ant  production  or  damping  caused  by  the  buoyancy  flux  at 
the  ocean  surface,  B0.  This  neglect  of  the  buoyancy  flux  is 
inconsistent  with  the  combined  contributions  to  the  mixed 
layer  buoyancy  that  will  be  associated  with  the  assumed  sur¬ 
face  heat  loss  of  500  W/m2,  the  assumed  evaporation  rate  of 
0.005  m/d,  and  a  predicted  freezing  rate  of  about  0.1  m/d.  As 
an  example  we  will  provide  a  one-dimensional  evaluation  of 
the  heat  budget  that  will  demonstrate  the  need  to  improve  (1). 

Including  B0,  the  recommended  form  for  the  entrainment 
parameterization  is 

w,  =  (ci«,3  ~  c2B0h1)/g*h1  (2) 

The  value  of  c,  is  about  2,  but  it  may  be  reduced  for  deeper 
mixed  layers  (for  comprehensive  reviews,  see  Zilitinkevich  et  at. 
[1979],  Garwood  [1979],  Gaspar  [1987],  and  Gallacher 
[1987],  In  tested  [ Martin ,  1985;  Gaspar,  1987]  parame- 
terizations  for  the  net  shear  production  of  TKE  less  dissi¬ 
pation  the  effective  value  of  c,  is  exponentially  reduced  on  a 
vertical  scale  of  ujf\ 

c,  ~  2  exp  (-fhjuj  (3) 

For  Hakkinen’s  numerical  simulation,  c,  would  be  reduced  to 
about  1.  The  value  of  c2  in  (2)  is  a  function  of  stability.  Stull 
[1976]  summarizes  values  of  c2  indicating  c2  ~  0.2  for  free 
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convection  (B0  <  0  and  ut  ~  0).  With  the  strong  wind  forcing 
in  the  Hakkinen  simulation,  c2  may  be  somewhat  larger. 

For  the  case  of  an  ocean  surface  just  beginning  to  freeze  (we 
will  assume  that  the  fractional  ice  cover  is  initially  small)  the 
surface  buoyancy  flux  has  three  components : 

S0  =  ~ocgQJ(pwcp)  -  flg{F  -  iWKS,  -  Sf)  -  pg(E  -  P)S , 

(4) 

Here  Qw  is  the  net  heat  loss  from  the  water  surface,  F  —  M  is 
the  freezing  rate  minus  the  melting  rate,  Sl  —  Sf  is  the  water 
surface  salinity  minus  the  ice  salinity,  and  E  —  P  is  the  evapo¬ 
ration  rate  minus  the  precipitation  rate.  To  contrast  the  differ¬ 
ences  between  (1)  and  (2),  we  shall  compute  B0  and  we  for 
Hakkinen’s  simulation  of  the  formation  of  new  ice  which 
occurs  in  a  region  that  is  nearly  horizontally  homogeneous 
during  the  first  day  [see  Hakkinen,  1987,  Figure  4c].  In  this 
region  we  evaluate  (2)— (4)  using  her  values:  E  —  P  =  0.005 
m/d,  F  —  M  =  0.1  m/d,  Sf  =  9.60  g/kg,  S,  -  34.60  g/kg,  S2  = 
34.95  g/kg,  Qw  =  500  W/m2,  =  75  m,  7]  =  - 1.89  K,  T2  =  0 

K,  a  =  0.025  x  10'3  K"1,  p  =  0.79,  /=  1.44  x  10"*  s_1, 
g  =  9.81  m/s2,  cp  =  4.186  x  103  J  kg'  1  K'  *,  and  pw  =  1028 
kg/m3.  We  find  that 

-CjBo/i!  =  0.2(2. 1  +  16.6  +  1.2)  x  10'4  m3/s3 

=  3.98  x  10'6  m3/s3  (5) 

for  the  sum  of  the  surface  heat  flux,  freezing,  and  evaporation 
contributions,  respectively.  Assuming  an  open  water  wind 
stress  of  3  dyn/cm2  (or  u„  =  0.017  m/s)  and  assuming  that 

<T  =  1, 

ci“*3  =  5.2  x  10“  6  m3/s3  (6) 

The  buoyancy  jump  at  the  base  of  the  mixed  layer  (reduced 
gravity)  is 

g *  =  Pg(S2  -  S,)  -  «g(T2  -  7])  =  2.29  x  10“ 3  m/s2  (7) 

Summing  the  buoyancy  flux  and  the  wind-stirring  contri¬ 
butions  to  the  entrainment  mixing,  (2)  is  used  to  compute  the 
entrainment  heat  flux  into  the  surface  layer: 

Q,  =  PwcpA7’w_  =  246  W/m2  +  188  W/m2  -  434  W/m2  (8) 

Conversely,  using  (1)  with  an  excessively  reduced  wind  mixing 
and  totally  neglecting  the  surface  buoyancy  flux,  we  calculate 
a  value  for  Qe  of  only  35  W/m2.  Clearly,  the  entrainment  rate, 
the  heat  and  salt  budgets,  and  the  rate  of  freezing  will  be 
significantly  affected  by  the  large  entrainment  heat  flux.  The 
freezing  rate  predicted  by  Hakkinen  of  about  0.1  m/d  was 
dependent  upon  a  net  mixed  layer  heat  loss  of  465  W/m2:  500 
W/m2  lost  to  the  atmosphere  less  35  W/m2  gained  by  en¬ 
trainment.  From  the  large  entrainment  heat  flux  predicted  in 
(8)  the  net  heat  lost  from  the  mixed  layer  should  be  reduced  by 
about  86%.  A  slightly  greater  entrainment  heat  flux  would 
prevent  freezing  entirely.  However,  an  important  feedback 
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Fig.  1.  A  positive/negative  thermodynamical  feedback  mechanism 
in  the  marginal  ice  zone. 


mechanism  cannot  be  neglected:  a  reduction  in  the  freezing 
rate  will  reduce  the  upward  buoyancy  flux  due  to  the  down¬ 
ward  surface  salinity  flux.  Thus  the  entrainment  rate  will  not 
be  quite  as  strong  as  that  computed  in  (5)— (8).  Nevertheless, 
we  recommend  that  (1)  be  replaced  by  (2)  in  Hakkinen’s  simu¬ 
lation  and  in  other  applications  of  the  complete  model. 

Buoyant  Damping  of  Mixing 
The  more  general  form  of  the  entrainment  equation  (2)  is 
especially  needed  when  precipitation  P  and/or  melting  M  are 
large  enough  to  cause  B0  to  be  positive  (downward  buoyancy 
flux).  Neglecting  this  effect  would  cause  excessive  deepening  to 
be  predicted.  This  may  explain  the  need  to  include  an  exces¬ 
sive  exponential  reduction  of  wind  mixing  in  (1).  Alternatively, 
if  the  recommended  form  (2)  for  the  entrainment  velocity  is 
used,  the  wind  forcing  may  be  balanced  by  the  buoyant  damp¬ 
ing,  and  the  mixed  layer  depth  will  asymptotically  approach  a 
steady  state:  the  Obukhov  length  scale,  i.e., 

hi  ~  ciu*3/cBo  (9) 

This  fundamental  process  in  mixed  layer  dynamics  is  not  pos¬ 
sible  with  (1). 


Feedback  Mechanisms 

The  alternative  entrainment  parameterization  suggested 
here,  equation  (2),  also  implies  a  positive/negative  thermody¬ 
namical  feedback  mechanism  in  the  MIZ,  as  shown  in  Figure 
1.  The  positive  feedback  mechanism  is  that  the  increased  (de¬ 
creased)  heat  loss  Qw  due  to  ice  melting  (freezing)  intensifies 
(weakens)  the  entrainment,  which  further  strengthens  (de¬ 
creases)  the  ice  melting  (freezing)  rate.  However,  a  negative 
feedback  mechanism  is  that  the  increased  upward  (downward) 
buoyancy  flux  caused  by  the  salinity  flux  due  to  ice  melting 
(freezing)  weakens  (intensifies)  the  entrainment,  which  further 
decreases  (strengthens)  the  ice  melting  (freezing)  rate. 

The  thermal  expansion  coefficient  of  seawater,  a,  is  very 
sensitive  to  temperature  and  is  small  at  the  freezing  point. 
Therefore  the  negative  feedback  mechanism  mentioned  here 
may  be  more  important,  in  general,  than  the  positive  feedback 
mechanism. 
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